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ABSTRACT 

We report the first detection of an extragalactic millimeter wavelength H2O 
maser at 183 GHz towards NGC 3079 using the Submillimeter Array (SMA), and 
a tentative submillimeter wave detection of the 439 GHz maser towards the same 
source using the James Clerk Maxwell Telescope (JCMT). These H2O transitions 
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are known to exhibit maser emission in star-forming regions and evolved stars. 
NGC 3079 is a well-studied nuclear H2O maser source at 22 GHz with a time- 
variable peak flux density in the range 3-12 Jy. The 183 GHz H2O maser 
emission, with peak flux density ~0.5 Jy {7a detection), also originates from the 
nuclear region of NGC 3079 and is spatially coincident with the dust continuum 
peak at 193 GHz (53 mJy integrated). Peak emission at both 183 and 439 GHz 
occurs in the same range of velocity as that covered by the 22 GHz spectrum. We 
estimate the gas to dust ratio of the nucleus of NGC 3079 to be 150, comparable 
to the Galactic value of 160. Discovery of maser emission in an active galactic 
nucleus beyond the long-known 22 GHz transition opens the possibility of future 
position-resolved radiative transfer modeling of accretion disks and outflows < 1 
pc from massive black holes. 

Subject headings: Masers — submillimcter — galaxies: active — galaxies: nuclei 
— galaxies: individual (NGC 3079) — techniques: high angular resolution 

1. Introduction 

Active galactic nuclei (AGN) that exhibit H2O maser emission at 22 GHz have been 
objects of particular interest since discovery that the maser emission can trace circumnuclear 
accretion disks, as in NGC 4258 (Miyoshi et al. 1995; Greenhill et al. 1995a,b). The presence 
of maser emission enables direct mapping of AGN dynamics within 1 pc of supermassive 
black holes via Very Long Baseline Interferometry (VLBI) (see reviews of Greenhill 2002; 
Maloney 2002). The additional discovery that H2O maser emission can be associated with 
shocked gas in jets and outflows has also provided a tracer for a second component in AGN 
(Gallimore et al. 1996; Claussen et al. 1998; Peck et al. 2003; Greenhill et al. 2003). 

Maser emission at 22 GHz is believed to be predominantly coUisionally-pumped at gas 
kinetic temperatures of Tj. = 400 - 1000 K, hydrogen densities of n(H2) = 10^ - 10^° cm~^, 
and dust temperatures T,^ < 100 K (e.g., Deguchi 1977; Elitzur et al. 1989; Yates et al. 
1997). In AGN accretion disks, gas heating is thought to arise from X-ray irradiation of the 
disk surfaces, whereas masers in jets and outflows may be heated via shock compressions 
(Neufeld et al. 1994; Kartje et al. 1999). 

Since the physical conditions that give rise to emission at 22 GHz also pump other H2O 

maser transitions, notably those at 183, 321 and 325 GHz (e.g., Deguchi 1977; Neufeld & 
Melnick 1991; Yates et al. 1997), 22 GHz maser emission is unlikely to be alone in AGN 
disks, jets and outflows. We note that, in galactic regions of star formation and evolved 
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stars, inversion of the 22 GHz transition is accompanied by inversion of other rotational 
transitions from the ground and vibrationally excited H2O states (e.g., Philhps et al. 1980; 
Waters et al. 1980; Menten & Melnick 1989; Menten et al. 1990a,b; Feldman et al. 1991; 
Menten & Melnick 1991; Melnick et al. 1993). In AGN, some of the masers should occur 
in broadly the same regions as at 22 GHz, but others could probe uncharted regions of 
the central engines of AGN, including regions at smaller radii. Line ratios of two or more 
maser transitions originating from the same gas would constrain radiative transfer models 
far better than is now possible. Previous searches for extragalactic millimeter/submillimeter 
H2O masers have been performed, but have been hampered by r.m.s. sensitivities worse 
than 1 Jy (Menten, private communication). 

NGC 3079 is an almost edge-on spiral Seyfert 2/LINER galaxy with optical heliocentric 
velocity of 1125 ± 6 km s~^ (e.g., Heckman 1980; Ho et al. 1997), at a distance of 16 Mpc 
(Sofue 1999). It hosts 22 GHz H2O maser emission at velocities between ~940 and 1350 
km s~^, which are believed to trace a disk-like structure at radii of 0.4 - 1.3 pc from a 
2 X 10® AIq central engine. The pc-scale maser disk is aligned with a dense molecular kpc- 
scale disk (Koda et al. 2002) and is located at the apex of a kpc-scale superbubble (Kondratko 
et al. 2005; Cecil et al. 2002). The bubble is believed to be inflated by a wide-angle outflow 
(Cecil et al. 2001) that is also responsible for the radio continuum lobes along the minor axis 
of the galaxy (e.g., Duric & Seaquist 1988; Trotter et al. 1998; Sawada-Satoh et al. 2000; 
Yamauchi et al. 2004). Within a parsec of the nucleus, the relativistic jet appears ahgned 
with the inner wall of the bubble and may be responsible for weak 22 GHz maser emission 
away from the disk (Kondratko et al. 2005, and references therein). All of these emission 
sources lie within the primary beam of the SMA observations reported here. 



2. Observations and Data Reduction 

2.1. SMA Observations at 183 GHz 

We observed the ground-state 3i3 220 transition of para-H20 (E„/A; ~200 K) at a rest 
frequency of 183.310 GHz (A: 1.6 mm) towards NGC 3079 using the SMA^ on 2005 March 
1. The compact configuration included six antennas, resulting in projected baselines of 8 - 
72 m. The SMA receivers operate in double-sideband mode with two 2 GHz sidebands that 
we tuned to 183.310 and 193.310 GHz (band-center). The correlator provided continuous 



^The SMA is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and the Academia 
Sinica. 
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coverage of each sideband (in 24 separate chunks of 256 channels) and 0.4 MHz channel 
spacing. The phase center of the observations was 0:2000 = 10'*01'"58*.53, (^2000 = 55°40'50''l 
(Cotton et al. 1999). The zenith opacity measured with the NRAO tipping radiometer 
located at the Caltech Submillimeter Observatory was r(225 GHz)=0.03 - 0.05 throughout. 

We cahbrated the data within the IDL superset MIR developed for the Owens Valley 

Radio Observatory (Scovillc 1993) and adapted for the SMA; the imaging was performed in 
MIRIAD. Data reduction took into account the effect of the pressure-broadened terrestrial 
H2O absorption line centered at 183.310 GHz. Fortunately, the rcdshift of the galaxy reduces 
the line frequency of the mascr transition by ~0.7 GHz, greatly lowering the corresponding 
atmospheric opacity. For bandpass calibration we used Jupiter and Saturn. Phase calibration 
was performed via frequent observations of the quasar 0923+392 about 18° from the phase 
center. We performed the amplitude calibration separately for upper and lower sideband 
data. For the upper sideband, we used 0923+392; for lower sideband data, we performed 
amplitude calibration on each of the 24 spectral chunks separately using Saturn and Jupiter. 
The calibration of the lower sideband both as a function of elevation and frequency was 
designed to reduce the effect of variation in terrestrial H2O absorption across the sideband. 
We derived the flux scale using observations of Callisto with an estimated accuracy of 20%. 
Single sideband system temperatures corrected to the top of the atmosphere were between 
400 and 850 K. 



2.2. JCMT Observations at 439 GHz 

We observed the ground-state 643 — > 659 transition of ortho-H20 (E„/A; ~1100 K) 
at a rest frequency of 439.151 GHz (A: 0.7 mm) towards NGC 3079 using the 15 meter 
JCMT on 2004 November 18 in r225=0.041 - 0.045 weather. The instantaneous bandwidth 
was 1840 MHz, corresponding to ~1200 km s~^, with a channel spacing of 1250 kHz or 
0.85 km s~^. The integration time on source was 2 hrs in beam-switching mode with a throw 
of 60". The average Tgys during the observations was 2650 K. We reduced the data using 
the Starlink spectral-line reduction software, SpecX. 

3. Results 

We detected millimeter (193 GHz) continuum emission with a peak signal to noise ratio 
(S/N) of ~10 (Fig. 1) at ctaooo = 10'*01"^57^80 ± 0.02, ^2000 = 55°40'46:'9 ± 0.3 (statistical 
uncertainties). Systematic uncertainties of (0f02,0'.'3) were estimated through measurement 
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of an apparent position for 1150+497 against the phase cahbrator 0923+392. These place 
the peak of the miUimeter continuum within la of the 22 GHz maser position. We partially 
resolved the nucleus in the north-south direction and the elongation is approximately parallel 
to the position angles of both the pc-scale accretion disk traced by 22 GHz maser emission 
and the kpc-scale molecular disk detected using e.g., CO (Koda et al. 2002). 

We detected 183 GHz H2O emission with a peak S/N of ~7 in both the real part of 
the amplitude spectrum and the interferometer phase spectrum (Figs. 1 and 2). The peak 
flux density is 0.55 Jy at 1017 km s^^. The line covers 9 frequency channels and overlaps 
a persistent Dopplcr component appearing in published 22 GHz maser spectra for dates 
1998 May 8, 2000 March 18 where it is especially prominent at ~2.5 Jy (both Hagiwara 
et al. 2002), 2001 March 23 (Kondratko et al. 2005) and 2002 April 12 (Braatz et al. 2003). 
We note possible detection of additional spectral features at 1208 (5 channels wide) and 
1265 km s~^ (4 channels wide). Both features appear in spectra of real amplitude and 
phase (near zero degrees, which indicates positions similar to the 1017 km s"-*^ emission), 
and both lie within the velocity interval of 22 GHz emission. Note that the correspondence 
of features between 22 and 183 GHz observations is not precise and is also not expected 
to be true for all features under certain pumping scenarios. The peak maser emission is 
at ^2000 = 10''0r'"57^75 ± 0.03, ^2000 = 55"40'46^'7 ± 0.4, which is within la of the 193 
GHz continuum peak. At our resolution, we cannot determine whether 183 GHz emission 
originates from the disk or outflow, or from both. We note that in addition to tracing a 
disk, a subset of 22 GHz maser emission is associated with the outflow within a few parsecs 
of the nucleus (Kondratko et al. 2005). 22 GHz emission at around 1000 and 1200 km s~^ 
is present both in the disk and outflow. Possible angular extension of the 183 GHz emission 
to the southeast is marginally detected (2a), which lies along the axis of the relativistic flow 
(e.g., Kondratko et al. 2005). We note that Hagiwara et al. (2004) detect OH absorption 
associated with the outflow at an overlapping velocity, and that we cannot rule out shocks 
in circumnuclear material as an origin for 183 GHz emission. 

We obtained a tentative detection {5a) of the 439 GHz maser at 1157 km s~^ (Fig. 3). 
With peak S^, ~2.5 Jy, 439 GHz emission is of comparable strength to typical values for peak 
22 GHz emission towards NGC 3079 and occurs within the 22 GHz velocity interval. 
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4. 



Discussion 



4.1. 



Millimeter Dust Emission 



The millimeter continuum of NGC 3079 is probably dominated by thermal emission 
from optically thin dust. As in the 1.2 mm observations of Braine et al. (1997) at 11" 
resolution, we find dust concentrated in the nuclear region of NGC 3079 (Fig. 1) and place 
an upper limit of 8 mJy {3a) on emission associated with the radio lobes ~1 kpc from the 
emission peak. Braine et al. (1997) estimates the dust temperature to be ~30 K at kpc radii. 
However, in two other AGN the presence of hotter dust at 100 times smaller radii, of ~300 
K, has been inferred from VLT observations (Jaffe et al. 2004; Prieto et al. 2004). Adopting 
a dust temperature of 100 K for NGC 3079 on intermediate scales and a dust opacity index 
/3 = 2, we estimate the enclosed dust mass of the inner 200 pc is M^ust ^ 10^-^='=°-^ Mq using 
the method described in Hildebrand (1983). Using CO Nobeyama observations by Koda 
et al. (2002), who find a gas mass of Mgas ~10^-^ Mq within the central 150 pc, we estimate 
a gas-to-dust ratio of 150, which is comparable to the Galactic value of 160 and to the ratios 
derived in bright local galaxies (Dunne et al. 2000). 



The 183 GHz emission we detect in the nucleus of NGC 3079 is likely from a maser 
process. It arises from one, and possibly more, narrow unresolved features in the same range 
of velocity as that of known 22 GHz maser emission. From radiative transfer models we 
know that the same high densities and temperature conditions that produce maser emission 
at 22 GHz also strongly invert the 183 GHz H2O transition. Noting that the 22 GHz 
feature at 1017 km s~^ varies between 0.4 - 2.6 Jy (Hagiwara et al. 2002), and assuming 
an ortho:para H2O ratio of 3:1, pumping models could reproduce likely line ratios (e.g., 
Yates et al. 1997). Narrow linewidths (<8 km s~^) at 183 GHz, small compared with the 
characteristic rotational and apparent random velocities for molecular gas in this nucleus 
(Kondratko et al. 2005), also indicate that emitting regions should be compact. Our beam 
size is too large to determine whether emission is thermal or maser based on brightness 
temperature arguments. For the emission to be thermal, with a brightness temperature of 
< 2000 K such that H2O molecules are not dissociated, the emitting region must be > 0''1 



4.2. 



183 GHz Emission 



(10 pc). 
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4.3. 439 GHz Emission 

439 GHz emission is probably also due to maser amplification. Emission occurs over a 
narrow velocity range (Fig. 3), and conditions for strong 439 GHz inversion are likely to exist 
in the galactic nucleus. Tiie 22 GHz spectrum peaks to the blue of the v^j^^ of 1125 km s~^, 
yet the 439 GHz spectrum peaks slightly above Vgys at 1157 km s~^ (Fig. 3). We ascribe 
this difference to the different pump mechanisms operating for each maser. Whereas 22 GHz 
emission is coUisionally-pumped and is increasingly quenched by local dust temperatures 
>100 K, the 439 GHz maser is predominantly radiatively-pumped and is most strongly 
inverted in the presence of = 300 K, = 400 K and n[H2) = 5 x 10^ cm~^ (Yates et al. 
1997). Dust temperatures of this order and higher are reasonable in the central parsecs of 
the nucleus (Jaffe et al. 2004; Prieto et al. 2004). Peak 439 GHz emission (1157 km s"^) 
is therefore likely to originate from a hot, dusty region around the central engine, possibly 
at radii <^1 pc, separated spatially from the gas emitting most strongly at 22 GHz (956 
km s"^). At 22 GHz, features at these velocities are indeed separated spatially by ~1 pc 
projected (Kondratko et al. 2005). The observations in Fig. 3 were taken on different dates, 
and we note that the 22 GHz feature near 1157 km s~^ is not present at all epochs (e.g., in 
Fig. 2). However, the feature line ratio in Fig. 3 of ~100 is accommodated by the vicinity 
of hot dust (Yates et al. 1997). 

We thank the all the SMA and JCMT staff. We thank S. Paine for opacity data, 
J. Braatz for 22 GHz data and G. Melnick for useful discussions. H. B. acknowledges 
support from the Emmy-Noether Program of the Deutsche Forschungsgemeinschaft (DFG- 
BE2578/1). 
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Fig. 1. — (filled color contours) Millimeter continuum image of NGC 3079 at 193 GHz 
observed with the SMA. Peak flux density is 27 mJy bm~^ with a la r.m.s. of 2.5 mJy bm~^; 
integrated continuum flux density is 53 mJy. (line contours) 183 GHz maser emission at 
1017 km s^^ (optical, heliocentric). Peak maser flux density in the image is 0.39 Jy bm~^. 
Positive emission is shown in solid lines and negative features in dashed contours. Contours 
are k x 2a, where k = —1, 1, 2, 3 and lcr=56 mJy bm~^. The black cross marks the position 
of peak 22 GHz H2O maser emission at 956 km s~^: asooo = 10'^01"^57^802 ± 0.001, ^2000 = 
55°40'47!'26±0.01 (Kondratko et al. 2005). The black solid arrow indicates the P.A. (~-10° 
east of north) of the pc-scale 22 GHz H2O maser disk (Kondratko et al. 2005) and of the 
kpc-scale CO disk (Koda et al. 2002). The dashed arrow marks the direction of relativistic 
outflow with uncertainty depicted by bracketing solid lines of P.A. 126° ± 5° (Trotter et al. 
1998; Sawada-Satoh et al. 2000; Middelberg et al. 2004; Kondratko et al. 2005). The arc to 
the east denotes the opening angle of the kpc-scale superbubble outflowing from the galactic 
nucleus (e.g., Cecil et al. 2002). The SMA synthesized beam has FWHM of 5"9 x 3'.' A (upper 
sideband) and 6''1 x 3''6 (lower sideband, shown top left), with a P.A. of 27° in both cases. 



- 11 - 




800 1000 1200 1400 1600 1800 



Velocity (Optical Heliocentric) (kms ) 

Fig. 2. — 183 GHz maser spectrum towards the nucleus of NGC 3079 corrected for atmo- 
spheric absorption obtained with the SMA. (top panel) Real part of the spectrum exhibiting 
a ~0.55 Jy feature at a velocity of 1017 km s^^. The channel spacing is 5 km s^^. Atmo- 
spheric opacity for a zenith angle of 50° is shown by the dotted line; peaks at 1850 and 2150 
km s~^ are due to ozone. Note that the noise rises at lower velocities due to the effect of 
atmospheric absorption, (upper bottom panel) 22 GHz spectrum obtained with a DSN 70 m 
antenna (Kondratko, private communication) between 800 - 1900 km s~^. (middle bottom) 
Real part of the 183 GHz SMA spectrum between 800 - 1900 km s~^. (lower bottom) Gor- 
responding fringe phase. Tentative maser features at 1208 and 1265 km s~^ are also noted 
with arrows. 
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Fig. 3. — (top) 22 GHz spectrum obtained using the GET from Braatz et al. (2003). (bottom) 
439 GHz H2O emission towards NGC 3079 observed with the JCMT. Peak flux density is 
2.5 Jy, and the la r.m.s. is 0.5 Jy. The channel spacing is 13.7 km s~^. 



